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Resumen 


Saucedo, H., Zavala, M., Fuentes, C. (enero-febrero, 2016). 
Estimación de parámetros de infiltración a partir de 
mediciones de avance de riego por melgas empleando las 
ecuaciones de Saint-Venant, y Green y Ampt. Tecnología y 
Ciencias del Agua, 7(1), 117-124. 


Se presenta un método para realizar la estimación de 
parámetros de infiltración con base en el empleo de las 
ecuaciones de Saint-Venant para describir el flujo del 
agua sobre el suelo, y la ecuación de Green y Ampt para 
representar el flujo del agua en el suelo. La estimación de 
los parámetros hidrodinámicos de conductividad hidráulica 
a saturación y presión en el frente de humedad, se realiza 
aplicando el método Levenberg-Marquardt. El modelo así 
obtenido, permite el ajuste de los parámetros hidrodinámicos 
a partir de datos de pruebas de avance de riego por melgas y 
de la textura del suelo. 


Palabras clave: ecuaciones de Saint-Venant, ecuación de 
Green y Ampt, método Levenberg-Marquardt. 


Introducción 


El objetivo del diseño de riego por melgas es 
conseguir la aplicación de una lámina de riego 
requerida por el cultivo de la manera más uni- 
forme posible, conservando una eficiencia de 
aplicación alta. El diseño del riego consiste en 
determinar el gasto de aporte y el tiempo durante 
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el cual se aplica dicho gasto en la cabecera de la 
melga para lograr la mayor uniformidad posible; 
es decir, en determinar el gasto óptimo para una 
longitud de melga específica y una caracteriza- 
ción hidrodinámica del suelo dada. Para realizar 
dicha caracterización hidrodinámica puede uti- 
lizarse, por ejemplo, información de una prueba 
de avance de riego, a partir de la cual se realiza 
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el procedimiento de ajuste u optimación de 
los parámetros, buscando minimizar el error 
obtenido entre las observaciones realizadas y las 
obtenidas mediante el modelo empleado para 
la simulación de la fase de avance del riego por 
gravedad. 

El objetivo del presente estudio fue desa- 
rrollar un modelo que permita hacer de forma 
automática la optimación de los parámetros de 
conductividad hidráulica a saturación del suelo 
y de la presión en el frente de humedecimiento, 
empleando una forma precisa de modelar el flu- 
jo del agua a superficie libre, para lo cual se hace 
uso de las ecuaciones de Saint-Venant; mientras 
que para la descripción del flujo del agua en el 
suelo, se utiliza la ecuación de Green y Ampt, 
la cual es una forma sencilla de modelar el fe- 
nómeno manteniendo bases físico-matemáticas 
en la representación. 


Flujo del agua sobre la superficie del 
suelo 


El flujo del agua con una superficie libre es 
modelado con las ecuaciones de Saint-Venant 
que resultan de la aplicación de las leyes de con- 
servación de masa y cantidad de movimiento. 
En una melga, la relación entre su ancho y el 
tirante de agua permite considerar las ecuacio- 
nes correspondientes al escurrimiento sobre una 
superficie de ancho infinito (Woolhiser, 1975). La 
ecuación de continuidad se escribe como: 


0 (1) 
ot 0x ot 


la ecuación de momentum se escribe en la forma 
recomendada por Saucedo, Fuentes y Zavala 
(2005): 


10q 2q0q q) gol 
a ls 22-00 
hat Ie Laly 1 ax EU 1)+B ot e 


donde q(x,t) =U(x,t)h(x,t) es el gasto por uni- 
dad de ancho de melga [L?T*]; x, la coordenada 
espacial en la dirección principal del movimien- 
to del agua en la melga [L]; t, el tiempo [T]; U, 
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la velocidad media; h el tirante de agua [L]; J, 
la pendiente topográfica de la melga [LL”); J, 
la pendiente de fricción [LL"]; V, = 01/0t, el 
flujo de infiltración [LT"], es decir, el volumen 
de agua infiltrado en la unidad de tiempo por 
unidad de ancho y por unidad de longitud de la 
melga; 1, la lámina infiltrada [L]; g, la aceleración 
gravitacional [LT”]; el parámetro adimensional 
$ = U,, /U, siendo a.=1-U,, /U, donde U,, es la 
proyección en la dirección del movimiento de la 
velocidad de salida de la masa de agua debido 
a la infiltración. 

La relación entre las variables hidráulicas q 
y h con la pendiente de fricción, denominada 
ley de resistencia hidráulica, es adoptada de 
acuerdo con Fuentes, De León, Saucedo, Par- 
lange y Antonino (2004), haciendo uso de una 
ley potencial de resistencia: 


] 6) 


h? 
qe A Dl 


donde v es el coeficiente de viscosidad cinemáti- 
ca del agua [L?T"*] y kes un factor adimensional; 
d es un parámetro adimensional que varía entre 
1/2<d<1 en función del tipo de flujo; los va- 
lores extremos d = 1/2 y d = 1 corresponden 
respectivamente al régimen de Chézy y al flujo 
laminar de Poiseuille. 

Las condiciones inicial y de frontera que 
deben sujetar a las ecuaciones de Saint-Venant 
para modelar la fase de avance son: 


q(x,0)=0 y  h(x,0)=0 (4) 
q(0,t)=4,, g(x,,t)=0 y h(x,,t)=0 (5) 


donde xt) es la posición del frente de onda 
para el tiempo t y q, es el gasto de aporte en la 
entrada de la melga. 

Para cerrar el sistema es necesario conocer la 
forma en que evoluciona en el tiempo la lámina 
infiltrada en toda posición sobre la melga, es 
decir, la ley de infiltración. 
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Flujo del agua en el suelo 


El modelo de Green y Ampt (1911) se establece 
a partir de la ecuación de continuidad y la ley 
de Darcy con las siguientes hipótesis: a) el perfil 
de humedad inicial en una columna de suelo 
es uniforme 0 = 6 ; b) la presión del agua en la 
superficie del suelo es hidrostática: y =h > 0, 
donde h es el tirante de agua; c) existe un frente 
de humedecimiento bien definido caracterizado 
por una presión negativa: y =-h,< 0, donde h, 
es la succión en el frente de humedecimiento; 
d) la región entre la superficie del suelo y el 
frente de humedecimiento está completamente 
saturada (flujo en pistón): 9 =8_ y K=K_ donde 
K, es la conductividad hidráulica a saturación, 
es decir, el valor de la conductividad hidráulica 
de la ley de Darcy correspondiente al contenido 
volumétrico de agua a saturación. La ecuación 
diferencial ordinaria resultante es la siguiente: 


dl 


h.+H)A9 
Ve=ak 1 4 Hp +A 
dt 


(6) 


I 


donde A6 = 6 — 6, es la capacidad de almacena- 
miento; Í es el volumen infiltrado acumulado 
por unidad de superficie de suelo o lámina 
infiltrada. 


Solución numérica 


Para simular la fase de avance del riego por 
melgas se hace uso del modelo desarrollado 
por Saucedo, Zavala y Fuentes (2011). La forma 
discreta de la ecuación de continuidad para la 
fase de avance se escribe como: 


o(h, +1,) 
¡+H(1- det 
(og, + 0) 3! (x, x;) +(1- 0), +1,) 
h, +1 
—[0q, +(1-0)q,, ]0t+(x, -x,,) o(h, +1,) | 
+1-0)(4,, +1,,) 


—[9h, +(1-q)h, +p1, +(1 1, ](x, -x,) 


+[qh, +(1-),, +91, +(1 PI, 6, -x))=0 (7) 


La ecuación de momentum guarda la misma 
forma discreta para las cuatro fases del riego: 


2qh(q, -q,)9t +(gh 7), -h,)0t 
+ lx, -x,)[09, +(1-0)q,-09,, -(1-0)9,] 
+ 30H (x,-x,) [0], +(1-0)], -J, ] (8) 


+porgh(x, -x,)[o1, +(1-0)1, 01, - (1-0), |=0 


En las ecuaciones (7) y (8), 0t es el paso de 
tiempo, wm y y son factores de peso en espacio y 
tiempo. El cálculo de los coeficientes se realiza 
con base en los valores pertenecientes al nivel 
de tiempo anterior 7= (1 — $), + d9, h=(1- 4) 
h + dh, (Saucedo et al., 2005). Se introducen las 
variables pequeñas (Strelkoff € Katopodes, 
1977): 6h, 6q,, 0h, y 0q,, tales que para las celdas 
interiores se tiene: h, = h, + 0h, h,=h, + 0h, q,= 
q,¡+ Óq, y q,=4,,+ 09, . Se sustituyen en las formas 
discretas de las ecuaciones de continuidad y mo- 


mentum y se obtiene un sistema de ecuaciones 
algebraicas, cuya solución permite avanzar en el 
tiempo los valores de las variables hidráulicas. 
El detalle de las entradas al sistema matricial 
puede consultarse en Saucedo et al. (2011). 

La ecuación de Green y Ampt (ecuación 
(6)) se resuelve numéricamente usando un 
método de diferencias finitas; el procedimiento 
se encuentra bien documentado en la literatura 
y puede consultarse, por ejemplo, en Burden y 
Faires (1985). 

Se ha utilizado un paso de tiempo constante 
At= 1.0 s para el acoplamiento de las ecuaciones 
de Saint-Venant, y Green y Ampt. La discretiza- 
ción utilizada para la solución de la forma com- 
pleta de las ecuaciones de Saint-Venant guarda 
semejanza con las reportadas en la literatura: 
Katopodes y Strelkoff (1977): At. =5s, Akanbi 
y Katopodes (1988): At_.__ =1s, Playán, Walker 
y Merkley (1994) At_. =2.12s. 


Optimización de parámetros 


Para estimar los parámetros hidrodinámicos K. 
y h, a partir de los datos registrados durante 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VII, núm. 1, enero-febrero de 2016, pp. 117-124 


gía y Ciencias del Agua, vol. VII, núm. 1, enero-febrero de 2016, pp. 117-124 


Tecnolo 


Saucedo et al., Estimación de parámetros de infiltración a partir de mediciones de avance de riego por melgas empleando... 


una prueba de avance, puede utilizarse el mé- 
todo Levenberg-Marquardt, en el cual el valor 
estimado de los parámetros en una iteración da- 
da se calcula teniendo en cuenta una expresión 
de la forma: 


TOY Tp) +M]ap =-J(Y K(p) (9) 


donde J(p) es la matriz Jacobiana relacionada 
con las variaciones de la función de posición del 
frente de avance x,=x/(t) respecto a cada uno de 
los parámetros a optimar (K, y h,); I es la matriz 
identidad; K(p), el vector de diferencias entre 
las posiciones del frente de avance observadas 
en campo y aquellas obtenidas con el método 
Levenberg-Marquardt en la iteración anterior; 
A se denomina parámetro de amortiguación y 
es importante en la convergencia del método, 
su valor puede determinarse, por ejemplo, con 
alguno de los procedimiento sugeridos por 
Griva, Nash y Sofer (2009). 
La matriz Jacobiana J(p) se calcula como: 


0Xp  0Xp 
0K, 0h, 
0X 5  0Xp 
0K., 0H,, 
1w= (10) 
0X 5  0Xp 
0K,, 0h, 


en esta ecuación, ¡=1, 2,..., n es el número de 
puntos que se registraron durante la prueba de 
avance, los cuales se recomienda sean diez a lo 
más, para reducir el esfuerzo de cómputo. La 
función de posición del frente de avance x,=x(t) 
es resultado de la aplicación del acoplamiento 
numérico interno de las ecuaciones de Saint- 
Venant, y Green y Ampt. La estimación de las 
derivadas se realiza numéricamente, a partir de 
la función x, = x(t) calculada para dos valores 
cercanos de los parámetros a optimar. 
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A manera de ejemplo, se toma el caso de una 
prueba de avance reportada en la literatura, 
registrada en el suelo franco de Montencillo 
por Fuentes (1992). Se dispone de la siguiente 
información: gasto unitario q, = 0.0032 m*/s/m; 
pendiente topográfica J, = 0.002 m/m, longitud 
de la melga L = 100 m, parámetros para la ley 
de resistencia d = 1 y k = 1/54; parámetro en 
la ecuación de cantidad de movimiento f = 2; 
valor inicial del contenido volumétrico de agua 
9, = 0.2749 cm*/crm*. El valor de la porosidad 
se asimila al contenido volumétrico de agua a 
saturación, de forma tal que 0, = 0.4865 cm?/ 
cm?. En la figura 1 se muestra la gráfica de los 
valores registrados del frente de avance. 

En las figuras 2 a 4 se muestra la aproxima- 
ción que se tiene para las iteraciones primera, 
tercera y séptima de la aplicación del método 
Levenberg-Marquardt. En la figura 5 se mues- 
tra la evolución del error cuadrático medio 
conforme avanzan las iteraciones del método. 
Los valores de los parámetros h,= 40.4 cm y 
K.=1.35 cm/h fueron obtenidos mediante la 
aplicación del procedimiento de optimización 
para reproducir datos de la prueba de riego 
efectuada en el experimento mencionado. 

Puede apreciarse que el número de iteracio- 
nes para lograr la convergencia es relativamente 
bajo, pero debe tenerse en cuenta que el mayor 
esfuerzo de cómputo se destina a la estimación 
de la matriz Jacobiana que aparece en las ecua- 
ciones (9) y (10). 


Conclusiones 


Se desarrolló un método para la estimación de 
parámetros de infiltración basado en el empleo 
de las ecuaciones de Saint-Venant para describir 
el flujo del agua sobre el suelo, y la ecuación de 
Green y Ampt para representar el flujo del agua 
en el suelo. La estimación de los parámetros 
hidrodinámicos de conductividad hidráulica a 
saturación y presión en el frente de humedad 
se realiza aplicando el método Levenberg- 
Marquardt. El modelo así obtenido permite 
el ajuste de los parámetros hidrodinámicos a 
partir de datos de pruebas de avance de riego 
por melgas y de la textura del suelo. 
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Figura 1. Valores observados en una prueba de avance reportada en la literatura. 
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Figura 2. Resultado de la primera iteración del método Levenberg-Marquardt, h,= 33.4 cm y K, = 1.58 cm/h. 
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Figura 3. Resultado de la tercera iteración del método Levenberg-Marquardt, h,= 38.3 cm y K, = 1.49 cm/h. 
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Figura 4. Resultado de la séptima iteración del método Levenberg-Marquardt, R?= 0.997, h,= 40.4 cm y K, = 1.35 cm/h. 
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Figura 5. Evolución del error cuadrático medio durante la aplicación del método Levenberg-Marquardt. 
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